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A spacecraft pushed by radiation has the major advantage that the power source is not included 
in the accelerated mass, making it the preferred technique for reaching relativistic speeds. There are 
two main technical challenges. First, to get signihcant acceleration, the sail must be both extremely 
light weight and capable of operating at high intensities of the incident beam and the resulting high 
temperatures. Second, the transmitter must emit high power beams through huge apertures, many 
kilometers in diameter, in order to focus radiation on the sail across the long distances needed to 
achieve high final speeds. Existing proposals for the sail use carbon or aluminum hlms. Aluminum 
in particular is limited by a low melting point, and both have low mechanical strength requiring 
either a distributed payload or complex rigging. Instead, we propose here a graphene sail, which 
offers high absorption per unit weight, high temperature operation, and the mechanical strength 
to support simple rigging to a lumped mass payload. For the transmitter, existing proposals use a 
compact high power source, and focus the energy with a large (hundreds to thousands of km) space- 
based lens. For optical drive proposals in particular, existing proposals require launch from the outer 
solar system, have severe pointing restrictions, and require difficult maneuvering of the beam source. 
Here instead we propose an active Fresnel lens operating at optical wavelengths, allowing smaller 
apertures of less mass, easier pointing with fewer restrictions, and probe launch from the inner solar 
system. The technologies for both the sail and the transmitter are already under development for 
other reasons. Worked examples, physically smaller and less massive than those suggested so far, 
range from a 1kg payload launched to 10% of the speed of light by a transmitter only 25 times the 
mass of ISS, to a larger system that can launch a 1000 kg payload to 50% of the speed of light. 


INTRODUCTION 

A spacecraft pushed by radiation has the major ad¬ 
vantage that the power source is not included in the ac¬ 
celerated mass. This makes it one of the very few tech¬ 
niques that might achieve the velocities needed for inter¬ 
stellar probes with human-scale travel times. However, 
the technical problems associated with both the sail and 
the source are daunting. The sail must be extremely light 
weight, yet either absorb or reflect photons well. Since 
high intensities are needed to achieve significant acceler¬ 
ation, it must operate at high temperatures. Finally it 
should be mechanically strong to simplify connecting it 
to the payload. 

The transmitter has a very different, but similarly dif¬ 
ficult, set of constraints. To focus radiation out to the re¬ 
quired distance, it needs to be physically huge while gen¬ 
erating coherent radiation across the aperture. It should 
be lightweight since it must operate in space and hence 
the components must be launched and maneuvered into 
position. It should allow a convenient place for launching 
probes, and allow transmitting in different directions, to 
launch probes to nearby stars and later send the energy 
needed for returning the data. 


PREVIOUS WORK 

A laser-driven aluminum sail was proposed by Forward 
in 1984[T]. The power source was a solar-powered laser. 


focussed by a 1000 km diameter Fresnel lens located at 15 
AU from the sun, and constructed of concentric rings of 
thin plastic film. The system performance was limited by 
the melting point of the aluminum sail. The type of laser 
source was to be decided later as technogy advanced, 
and the details of aiming the laser beam and keeping it 
focussed on the probe were not addressed. 

A lightweight, microwave-driven probe was introduced 
as Starwisp[2], with an aluminum mesh driven by mi¬ 
crowave radiation. The microwaves were to be generated 
by a solar power satellite, and used only incidentally for 
accelerating spacecraft. The wavelength was not optimal 
for this application and required a truly enormous lens, 
much larger than the Earth. Furthermore, the payload 
was tiny, only 4 grams, and had to be distributed across 
the sail in even smaller units, a technology not yet mas¬ 
tered and difficult even at a conceptual level. Finally, 
later analysis, first informal [3] and then formal[l], de¬ 
termined the aluminum mesh would melt under the pro¬ 
posed illumination. Scaling the power density down to 
the point the mesh was not melted reduced the available 
acceleration dramatically. 

All of these issues were addressed by Landis [4]. After 
recognizing that practical materials will operate in ab¬ 
sorption and not reflection, he advocated a new figure of 
merit for potential sail materials. For a sail dominated 
by absorption, the thrust is P/c, where P is the power 
absorbed. This energy must be lost by radiation, which 
grows as T^. Therefore the highest achievable accelera¬ 
tion scales as T^/m, where T is the operating tempera- 
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Work 

Starwisp 

microwave 

Landis 99 
microwave 

Forward 84 
optical 

Landis 99 
optical 

Power 

10 GW 

56 GW 

65 GW 

448 MW 

Lens size 

50,000 km 

125 km 

1000 km 

2600 km 

Wavelength 

3 cm 

3 mm 

1 um 

400 nm 

Sail size 

1 km 

100 m square 

3.6 km 

4.08 m 

Payload 

4 gram 

80 gram 

333 kg 

6 gram 

Acceleration 

115 g 

2.5 g 

0.036 g 

4.17 g 

Sail temp 

melts 

2333° K 

600° K 

1583° K 

Acc. Time 

10 hours 

18.5 hours 

3 years 

8.5 days 

V at cutoff 

0.1 c 

0.005 c 

0.11 c 

0.1 c 


TABLE I: Summary of previous work 


ture and m the mass per unit area, and the best perfor¬ 
mance is obtained with lightweight materials with high 
melting points. Carbon films appeared the best choice in 
these studies. 

Next, Landis proposed the use of shorter wavelength 
microwave illumination (3mm instead of 3cm), which re¬ 
duced the lens size to 125 km. Finally, he used a larger, 
more practical, lumped payload of 80 grams. This re¬ 
quired a complex fractal design for the rigging to connect 
it to the sail without overstressing any portion of the thin 
membrane. 

Landis also briefly mentioned an optical sail made of 
sapphire. The sail is 57 nm thick to make it A/4 at 
400 nm incident light, to optimize reflectance. It reflects 
26% of the light that falls on it while absorbing very lit¬ 
tle - the absorption was too low to for him to measure 
experimentally, and he assumed 0.5%. The system was 
only partially fleshed out, but has several numerical in¬ 
consistencies. From the stated parameters (laser = 448 
MW, reflection = 26%, acceleration = 43.4 g), the to¬ 
tal mass of the sail -|- payload is 1.8 grams, which seems 
very small, particularly since there is a discussion earlier 
in the paper of whether a payload can be reduced to 5 
grams. The payload mass is not explicitly stated, but 
the text states the total mass is 2.3 times the sail mass, 
citing Forward’s 1984 paper. However, that paper used a 
ratio of 3 (“roughly one-third each of sail, structure, and 
payload”). Finally, the time of acceleration was 8.5 days, 
and the final velocity was 10% c, but that combination 
gives 4.17 g, not 43. Assuming that at least a Starwisp 
size payload is needed, we will assume what was meant 
was 4.17 g, 8.5 days of acceleration, 18 grams total mass, 
and 6 grams payload to 0.1c. At the stated power den¬ 
sity of 34 MW/m^ this gives a sail 4.08 m in diameter. 
Although the required power is indeed low, it requires a 
huge lens (2600 km in diameter) to focus on such a small 
sail out to the required 75 AU. 

An overview of these proposals is shown in Table |lj 

Graphene sails were mentioned briefly by Benford in 
[2 , where their high temperature capabilities were noted. 
This work also tried to minimize the cost of the transmit¬ 



FIG. 1: Absorption as function of number of layers, from 
[9]. Absorption in the visible is on the right edge (and futher 
right) of this figure. Absorption is proportional to the number 
of layers, at least up to 8, and shows no sign of deviating from 
proportionality. 


ter, with different techniques than pursued here. Matloff 
analyzed the use of graphene in solar sails [^[7], noting 
the same overall advantages as mentioned here, though 
different criteria are needed for this application. Bible et 
al.[S] also analyze a large scale phased array laser trans¬ 
mitter, with many fewer but much larger elements than 
used here. This work assumes a reflective, not absorbing, 
sail. 

OVERVIEW 

This paper continues these trends to create a still more 
“practical” solution. The first innovation is the use of 
graphene for the sail. This provides three advantages 
over a carbon film. The absorption is very high per unit 
thickness (in the optical, at least), it’s mechanically very 
strong, which simplifies the rigging, and it’s capable of 
high temperature operation. The second innovation is 
the use of a phased array, rather than a lens, for the op¬ 
tical drive. This offers flexibility in aiming, probe launch 
from near-Earth, and does not require moving the source 
during operation. 

Graphene consists of one or more layers, each a two- 
dimensional hexagonal array of carbon atoms. Since 
the studies of Landis, it has achieve considerable promi¬ 
nence. Its invention garnered the Nobel prize in physics 
in 2010 HH], and the award included a good summary 
of its properties. Single layer graphene has a den¬ 
sity of 0.77mg/m^, and an optical absorbance of 2.3%, 
independent of wavelength (in the visible spectrum). 
Graphene is 0.345 nm thick, and has physical strength 
of 42 N/m, well above that of other materials of similar 
thickness. The exact melting point is hard to determine 
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experimentally [TT], but appears to be about 4900° K[T^. 

Consider a sail made of single layer graphene with a 
working temperature of 4000° K. It will radiate mostly 
in the optical, so its emissivity will be the same as the 
absorbance, or 2.3%. Thus it radiates (2 sided) about 
2aeT^ or about 668 Kw per square meter. If it absorbs 
this all from one side, it applies a force of 668 kw/c/m^, 
or 2.2 • 10“^ newtons/m^. Acting on the mass of 7.7 • 
10“^ kg/m^, this gives an acceleration of 2890 m/sec^, 
or about 300 g. Naturally the addition of a payload, 
and the limits of available power, will reduce this figure 
considerably, but it shows the potential. 

Furthermore, a graphene sail is strong enough to use 
simple rigging to propel a large, dense object, which is 
not possible with Starwisp. Starwisp was constructed 
of O.I micron wires on .3 mm spacing. Assuming high 
strength aluminum of 440 MPascal strength, this gives a 
strength of 0.0012 n/m. Under the proposed illumination 
of 8.6 kw/m^, perimeter rigging (like a parachute) would 
exceed the tensile strength of the sail at a radius of 100 
m, with a total force of 0.7 newtons. Graphene, at over 
10,000 times the tensile strength of the aluminum grid, 
allows sails of almost any size to be rigged only at the 
rim. 

Graphene comes in both single and multi-layer forms, 
and the use of multi-layer graphene helps considerably. 
Each additional layer absorbs 2.3% of the remaining light, 
so for small number of layers the absorption is nearly 
proportional to the number of layers [13] ■ This has been 
measured directly for up to 8 layer graphene (see Fig[Tl 
from i), where the property holds without any sign of 
converting to bulk graphite, in the near-infrared (1 ev) 
and shorter visible wavelengths 

As additional layers of graphene are added, the mass 
goes up, but so does the absorption (and hence thrust 
from constant flux). The temperature is unchanged since 
emissivity rises at the same rate as absorption. For a 
graphene sail alone, a single layer is optimal, since adding 
another layer doubles the mass, but does not quite double 
the absorption. But if the payload is a significant fraction 
of the total, then the overall performance goes up as the 
number of layers increase. 

The limit to this process is when each additional layer 
adds less absorption than it does mass, or when the prop¬ 
erties of multi-layer graphene change to those of bulk 
graphite[3|. Assuming the trend in Fig continues to 
many layers, and the sail mass is equal to the payload (a 
reasonable approximation, as shown later), then the op¬ 
timum number of layers is 54, for a thickness of 18.6 nm 
and absorption of 72%. At this thickness an additional 
layer contributes 0.92% to both absorption and mass. 


(a) 

Material 

Aluminum 

Sapphire 

Graphene 

(b) 

Thickness 

16 nm 

57 nm 

17 nm 

(c) 

Mass/m^ 

43 mg 

226 mg 

42 mg 

(d) 

Reflectance 

82% 

26% 

0% 

(e) 

Absorbance 

13.5% 

1.7- 10“’’ 

72% 

(e) 

Thrust/W, 

1.77/c 

0.52/c 

0.72/c 

(f) 

Emissivity 

0.06 

0.10 

0.72 

(g) 

Melting pt 

933° K 

2303° K 

4900° K 

(h) 

Area for 333 kg 

7.7 km^ 

1.4 km^ 

7.9 km® 

(i) 

Radius for 333 kg 

1570 m 

685 m 

1589 m 

(j) 

Can radiate/m^ 

1.02 KW 

63 KW 

9.3 MW 

(k) 

Incident power/m^ 

7.54 KW 

136 MW 

12.9 MW 

(1) 

Power in beam 

58.4 GW 

201 TW 

102 TW 

(m) 

Thrust 

346 N 

348 KN 

246 KN 

(n) 

Total mass 

1000 kg 

1000 kg 

700 kg 

(o) 

Acceleration 

0.35 m/s^ 

348 m/s^ 

351 m/s® 

(P) 

Time to 0.1c 

2.73 yrs 

1.00 day 

0.99 day 

(q) 

Distance to 0.1c 

8690 AU 

8.65 AU 

8.56 AU 

(r) 

Total energy (J) 

5.07- 10^® 

1.73- 10^® 

8.74- 10^® 

(s) 

Transmitter radius 

253 km 

576 m 

246 m 


TABLE II: Comparison of sail material options. Data for the 
alnminum sail from [I]. Data on the sapphire sail from [4], 
except for transparency which is estimated from data in m- 
Below the line are compnted properties for the case of a 333 
kg payload to 0.1c . 


Comparison of sail materials 

The three materials proposed for optical sails are alu¬ 
minum, sapphire, and graphene. Their properties are 
compared in Table [TT] One way to compare them is to 
examine the results of using each material to accelerate 
a canonical 333 kg payload to 0.1c, an example already 
used by Forward. We will start by assuming an optimal 
transmitter for each sail. 

The calculation proceeds as follows. From both For¬ 
ward and Landis, and later in this work, the sail mass 
should roughly equal the payload mass. So the first step 
is to find the size of a 333 kg sail of each material, entries 
(h) and (i) in the table. Next, using the emissivity and as¬ 
suming 2-sided radiation, figure out how much power the 
sail can radiate, per m^, without exceeding 2/3 of melt¬ 
ing temperature. This is shown in line (j). The difference 
is dramatic, driven by the 1770-fold more power that can 
be radiated at the higher temperature with the higher 
emissivity. This difference in absorbed power dominates 
all the aluminum-graphene differences that follow. 

The next step is to use the absorption to figure out the 
incident power the sail can support. For sapphire, this 
gives an enormous value due to its very low absorption. 
However using the full value tolerable by the sail would 
doubtless overheat the payload. Even for the graphene 
sail, the illumination is already 10000 suns, requiring a 








4 


roughly 99.99% reflective payload to avoid overheating. 
We will give the sapphire sail the benefit of the doubt, 
and assume that a 99.999% reflective mirror over the re¬ 
quired 10% bandwidth is possible.^ This limits the beam 
intensity to about 100,000 suns, or 1.36 • 10® W/m^. Re¬ 
sults are shown in line (k). 

Given the incident intensity, the area, and the efh- 
ciency, the thrust can be computed, line (m). This acts 
on the mass, and here the lesser strength of aluminum 
and sapphire come into account. We assume the mass 
of the supporting structure is equal to the payload mass 
for aluminum (as assumed by Forward and Landis). We 
also assume this for sapphire, though this again might 
give it the benefit of the doubt since the huge intensity 
in that case would likely require sapphire rigging. And 
for graphene, its strength allows simpler rigging (we will 
show this later), so we assign it 10% of the mass of the 
sail. The total mass is shown in line (n). 

At this point we have the acceleration (o), from which 
we compute the time (p) and distance (q) needed to reach 
0.1c. From the distance, and the size of the sail, we derive 
the radius of the transmitter needed (s) (assuming A = 
500 nm in all cases, though this too is an approximation). 
From the time of acceleration, and the beam power, the 
total energy expended (r) can be computed. 

As mentioned earlier, the dominant factor in these cal¬ 
culations in the melting point, since it is raised to the 4th 
power. The aluminum sail has the best efficiency - it gen¬ 
erates more thrust per photon and has comparable mass 
to graphene (and much less than sapphire). However, it 
has a low maximum temperature, limiting the peak ac¬ 
celeration, hence requiring a long distance, a long time 
(almost three years), and a proportionally large lens - 250 
km in radius in this case. (This is half the radius found 
in Forward since we are using half the wavelength, but 
otherwise the calculations largely agree.) In return for 
the slow acceleration, the source power is much smaller, 
“only” 58 GW. 

The sapphire sail is an intermediate case. It has 
the highest total energy needed, and the highest power 
source, due to its low efficiency. But the acceleration is 
high due to the enormous intensity it can support, so the 
acceleration distance is short, comparable to graphene. 
But the lens is bigger, twice the size of graphene, since the 
higher mass per unit area means a smaller sail. Finally, 
two other factors would need to be taken into account 
for a sapphire sail. The reflectivity is tuned to light of a 


^ The 10% bandwidth requirement is caused by redshift as the sail 
accelerates. In theory, this could also be dealt with by changing 
the laser wavelength as the probe accelerates, to keep the wave¬ 
length constant at the probe. Then the mirror needs only be 
99.999% reflective at one wavelength, which has been done. But 
a super high power, efficient, and tunable laser seems at least as 
much of a stretch as a high bandwidth mirror. 


(b) Active 

(a) Reference 



radius sail 

FIG. 2: Schematic diagram of system. 


specific wavelength arriving perpendicular to the surface, 
and will work less well as redshift increases or if the sail 
is not at right angles to the incoming radiation. 

The graphene sail’s combination of strong absorption, 
high temperature operation, and light weight yields the 
smallest lens, as small as 246 m in radius. (In practice, 
other considerations such as energy density may impose 
larger minimums) The source power is high, 102 TW, 
but allows a probe launch once per day, rather than once 
every three years. 

The possibility of a smaller lens allows the considera¬ 
tion of more complex lens architectures. Since the lens 
used by Forward was large, it was essential to keep the 
mass per unit area as low as possible, and hence it was 
a passive structure only a half wavelength thick on the 
average. But if the lens size can be reduced, perhaps by 
an order of magnitude, then a much thicker and more 
complex lens can be considered, and still have less mass 
than the large thin lens. 

As a side note, if you want only to accelerate the sail, 
and are not concerned about the payload, then sapphire 
offers some incredible performance. The absorption of 
visible light by sapphire is extremely low. Even at high 
temperatures (1500° K) the absorption is only about 0.03 
per cm. m So a 57 nm film should only absorb about 
1.7-10“^ of the incident light. This would allow an enor¬ 
mous incident intensity of 411 GW/m^. Reflected with 
26% efficiency, this generates a force of about 712 N/m^, 
which acting on a mass of 226 mg/m^, creates an accel¬ 
eration of 2.68 • 10® m/s^, or about 273,000 g. If this 
acceleration could be sustained, it would reach the speed 
of light in less than two minutes. Unfortunately, using 
this theoretical performance for anything more than ac¬ 
celerating the sail alone seems impractical. Any payload 
attached to the sail will be instantly vaporized unless it 
is either as transparent as the sapphire or far more re¬ 
flective than any known material. 
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OVERALL ARCHITECTURE 

Existing proposals require two large space-based struc¬ 
tures - beam generation and the lens. For example, the 
proposal of Forward starts with a 65 GW solar powered 
laser. At 25% efficiency and in Earth’s orbit, this in¬ 
volves 191 km^ of collecting area. This generates the 
beam which is then focussed by the 1000 km diameter 
lens in the outer solar system. 

Given the smaller lenses that are possible with 
graphene based sails, it is perhaps possible to combine 
these two functions into one structure, one that receives 
incoherent raw sunlight on one side and emits a coherent 
optical beam on the other, based on an optical reference 
signal. This is in essence a Fresnel lens, using a thin 
flat plate with local modification of phase to achieve lens 
operation. It is active since it adds power to the refer¬ 
ence signal, and since the local phase changes are un¬ 
der computer control. This allows relaxing the physical 
constraints on the lens, varying the focal length during 
operation, and control of the beam profile. 

Such a system might be built as shown in Fig A 
phase reference laser at (a) illuminates the back of the a 
multi-km radius transmitter, (b). This transmitter uses 
sunlight (thick arrows), solar cells, and phase controllable 
semi-conductor lasers to create a multi-GW coherent op¬ 
tical beam at optical wavelengths. This is partially ab¬ 
sorbed be the graphene sail at (c), which pulls the lumped 
payload (e) through the rigging (d). Since the rigging is 
at the edge, but the thrust is distributed, the sail be¬ 
comes convex which is needed for passive centering of 
the sail on the transmitted beam[T3. The acceleration 
of the sail and payload continues until the transmitter 
can no longer focus a significant fraction of its power on 
the sail. 

This approach solves several problems inherent to the 
configuration with a separate power source and passive 
lens. First, to keep the mechanical tolerances of the 
lens reasonable, it was built with an extremely long focal 
length - many AU for a 1000 km diameter lens. In the 
paper, the lens was proposed to be located 15 AU from 
the sun. At this distance, the lens-transmitter combina¬ 
tion can only cover a very small part of the sky, within 
a few degrees of the sun-lens line. For the specified 15 
AU distance, this line will orbit the sun each 58 years. 
This could require additional lenses should it be neces¬ 
sary to send a similar beam to the same target star at a 
particular time many years later, perhaps to power data 
return. 

The next problem is that the probe(s) must be 
launched from the outer solar system, presumably de¬ 
livered there by some more conventional means. This is 
closely related to the third problem, that the laser source 
needs to be located at one specific location in space. As 
with any conventional lens, as the probe moves away, the 


needed location for the source moves closer to the lens. 
This results in a tradeoff, where a larger focal length lens 
gives less required movement of the source, but the clos¬ 
est attainable focus is further away. 

For example, the design might want the laser to have to 
move at most 1 AU as the probe accelerates, as might be 
needed if the laser is solar powered and must remain close 
to the sun. Assuming a 15 AU lens location, this requires 
a 14 AU lens focal length which can only focus as close 
as 210 AU from the lens, or 225 AU from the sun. This 
is much further out than any man-made object to date 
(Voyager is at 131 AU as of 2015, after almost 40 years of 
flight.) Even with this large starting distance, the light 
source must accelerate up to 7 km/sec, then decelerate, 
during the launch of the probe. A shorter focal length 
lens can offer the probe a starting spot closer to the sun, 
but requires more movement and higher speeds from the 
laser. If you want to start at only 30 AU from the sun, 
for example, the focal length needed is 7.5 AU, and the 
focal point moves more than 7 AU (making solar power 
difficult) at a peak speed of 270 km/sec (well above the 
Sun’s escape velocity). 

An active lens solves all these problems. Since the 
active lens focal length is under computer control and can 
be changed at will, the lens could be built at a Venus orbit 
(for example) and the probes constructed and launched 
near Earth. The lens can cover a wide portion of the 
sky - up to a 45 degree offset from the sun-lens line can 
be accomodated with a degradation in power. In the 
worst case, 8 such lenses, four in the orbit of the ecliptic 
and four in sun polar orbit, could cover the whole sky at 
any time. 

The following sections work through some of the 
needed calculations, then size and performance estimates 
are given for a variety of possible interstellar missions. 

FOCUSSING DISTANCE 

How far away can a transmitter of radius K focus a 
beam onto a sail of radius r, using radiation of wavelength 
A? As in [T] , we define this distance as that where the sail 
just fills the main lobe of the diffraction limited circular 
beam. At this point 84% of the total beam energy falls 
upon the sail. If the transmitter is radius i?, the sail 
radius r, then this happens at distance 

4ri? 

^ 2.44A 

This distance is known as “cutoff” in the literature. 

After cutoff, the final speed can still be increased, but 
the acceleration now falls as 1/s^, where s is the distance. 
In the classical case, if continued forever, the energy will 
eventually double and the speed increased by a further 
•\/2. This result is stated, but not derived, in [5]. Perhaps 
the simplest derivation is via calculating the work done 
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on the craft by the photons. Assume we have force Fq 
applied over a distance S'o, with force Fo(S'o/s)^ there¬ 
after. Then the energy imparted up to cutoff is simply 
FqSq, and the work up distance s, greater than cutoff, is 

E{s) = FoSo+ f FoiSo/sfds 

JSo 

Integrate to get 

E{s) = FoSo{2 - So/s) 


s(t)=ct-^(l-e^‘) (3) 

Since the distance D to which the beam can focus is un¬ 
affected by relativistic corrections, to find the final rela¬ 
tivistic speed we invert (3), setting s{t) = D, to find the 
time t of last drive. We plug this time t into equation (2) 
to get the final speed. 

DISK STIFFNESS 


and the energy is at most doubled by acceleration after 
cutoff. However, acceleration after cutoff is an inefficient 
use of a very expensive beam generator, and hence for 
the purpose of this paper we assume drive stops when 
the beam can no longer focus down to the size of the 
sail. 

It should be noted that there is confusion on this point 
in the literature. Although Forward correctly states (see 
[2], equation (12)) that it is energy that can be doubled 
by post-cutoff acceleration, many of the examples in [5], 
[1] , and [T] assume instead it is the velocity that is dou¬ 
bled. This results in some unjustifiably optimistic final 
velocities. 


RELATIVISTIC CALCULATIONS 


Since the whole point of this operation is to attain 
relativistic speeds, relativistic calculations are necessary. 
Consider a mass Mg propelled by radiant power Po, 
which is entirely absorbed. This induces a force Fg = 
Pg/c and and initial acceleration Ag = Now as the 
mass picks up speed, the acceleration goes down in two 
ways. First the incoming photons are red-shifted, so the 
force goes down 


E{v) 


^ j l - v/c 
c y xfv/c 


In addition, the mass goes up 


M{v) 


Mq 

\Jl — jc? 


and hence the acceleration (in the rest frame) is 


A{v) 


MqC 


y 1 -I- v/c 


Pq 

MqC 


(1 - v/c) 


This sets up a differential equation, and assuming a(0) = 
?;(0) = s(0) = 0, the solutions are: 

ait) = Age^‘ (1) 


To keep the sail in the beam without active control, it 
needs to be convex as seen by the beam[T5]. This could 
be achieved by attaching the rigging to the rim of the 
disc, then rotating the disk to keep it from collapsing. 

A graphene disk can rotate at high speeds. A uniform 
disk of density p, radius R, and Poisson’s ratio z/, rotating 
with angular speed w has a maximum stress cr |I6j 


O’M = (3 -I- v) 


puj^R^ 

8 


Working backwards from the properties of graphene (tr = 
10^^ Pascal, density = 1300 kg/m^, Poisson’s ratio of 
0.25), we find (for example) that a disk of 1km radius can 
rotate with a speed of 16 radians/sec. This gives a speed 
at the edge of 16,000 m/s, and an acceleration of 256,000 
m/s^, or about 25,000 gs. This is much larger than the 
linear acceleration and keeps the disk from collapsing. 
The hoop stress at the outer edge is given by: 

puj^R^ 

an = {1- v )—— 

or, comparing to the maximum stress above 

(1 - zz) 8 


o’/f = ctm 


(3 -k z/) 4 


Using the above values, the maximum hoop stress is 
therefore about 4.6 • 10^^ Pascal. 

Assuming the rigging is at a 45° angle, the sail will 
experience a force trying to collapse the sail that is equal 
to the thrust. The outward hoop stress from the spin¬ 
ning sail needs to exceed this force to prevent the sail 
from collapsing. The calculations below Hnd the rim size, 
the distance from the edge where the integrated outward 
stress exceeds the inward stress from the rigging, assum¬ 
ing the disc is rotating as fast as possible. In all cases this 
is a very small fraction, much less than 1%, of the sail 
radius. In practice, the sail would not rotate at its maxi¬ 
mum possible rate, but even so the rigging only needs to 
attach to a small area near the rim, and the rest of the 
sail will be in pure tension. 


OPTIMUM SAIL SIZE 


u(t) = c(l — e “”*) 


This section demonstrates, that at least in the classical 
(2) case, the optimum final velocity (for constant thickness) 
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is achieved when the sail has the same mass as the pay- 
load. Changing the sail size has two effects - it changes 
the total accelerated mass, and it changes the distance 
to which the transmitter can focus. It does not change 
the thrust since the entire beam is incident upon the sail 
independent of the sail size. So we have 

v{r) = y/2a{r)s{r) 

If we set P to be the payload mass, E the absorbed power, 
and Mq the mass per unit area times tt, then we get 


, , L E 4rR 

“ Y ^ ■ c{P + Mor2) ■ 2.44A 

Now we need to find the optimum r. First, since the 
square root of a positive number is a monotonic function, 
we can just find the optimum of 

^ E ArR 

^ ' c[P P Mor2) ■ 2.44A 

next, £1, i?, c, and A are just multiplicitive constants, and 
do not change the optimum r. So we just need to find 
the r that maximizes 

r 

P P Mor2 

The derivative with respect to r is 

(P-f Mor^) • 1-r •2Mor 
(P-f Mor2)2 

This can only be 0 if the numerator is 0, so we set 
P -P Mor^ - 2Mor^ = 0 

This is 0 when = P, or in other words when the 

mass of the sail equals the mass of the payload. From 
the numerical evidence, the same appears true in the rel¬ 
ativistic case, where v = y/2as is replaced by inverting 
equation 3 and then inserting the time into equation 2. 
The analytic proof of whether this is precisely true is left 
to the interested reader. 

THE ACTIVE LENS 

How can an active lens such as proposed here be con¬ 
structed and aimed? Unlike the sail, the active lens had 
many possibilites for its construction. The arguments 
here are therefore to show that it is plausible that it could 
be built, not a blueprint for construction. The main idea 
is to combine a high efficiency solar cell [17] with a high 
efficiency laser[T8]. To minimize the size, mass, and cost 
of the lens, the power output per unit area should be as 
high as practical. This implies high illumination on the 
solar cells, either by locating the lens closer to the sun 



FIG. 3: Detail view of active lens. Module C computes the 
required phase shift based on the measured location of the 
pixel and the desired beam direction. Module 0 phase shifts 
the reference beam, which is then amplified by (or locked to) 
the VCSEL. The VCSEL gets its power from the solar cell on 
the sunward side. The distance to the lens is not to scale, as 
the exit divergence of VCSELs is 5-10°. 


than the Earth, or perhaps using large thin reflectors to 
concentrate additional light onto the array. The limiting 
factor is when the temperature of the array becomes too 
high for reliable operation. We pick a reference insolation 
of 4 kw/m^. Since good solar cells convert >40% of light 
into electricity, and good lasers have more than >60% 
efficiency, this gives an overall efficiency of perhaps 25%, 
generating Ikw/m^ of drive power. If we assume that of 
the 4kw/m^ incident, Ikw/m^ is reflected and Ikw/m^ is 
the output beam, then this leaves 2 kw/m^ which must 
be re-radiated as waste heat. If the lens is also an effi¬ 
cient radiator in the IR, this leads to a temperature of 
91° C. 

The array consists of pixels spaced a few hundred wave¬ 
lengths apart. Each pixel continually adjusts its phase to 
generate a coherent output beam, based upon its physical 
location and the desired output direction. The reference 
beam arrives from the rear. It is phase shifted on a per- 
pixel basis, then amplified and emitted. Since the pixels 
are many wavelengths apart (perhaps 200 wavelengths, 
for a 100 micron spacing and 500 nm light), each pixel 
must emit only over a narrow angle (roughly 1 over the 
spacing in wavelengths) to avoid serious grating losses. 
This means the disk must be physically aimed close to 
the required direction. 

Multi-junction III-V compound cells are currently fa¬ 
vored for their high efficiency, and small thickness com¬ 
pared to silicon cells [TH). Assuming that the lens is dom¬ 
inated by the mass of the 10 micron thick solar cells with 





































Vertex 

spacecraft 
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FIG. 4: Aiming the lens - the big picture. The four filled 
squares fly in formation, defining a reference tetrahedron by 
keeping constant length on the six thick dashed lines. Each 
element of the lens computes its coordinates by comparing the 
modulation and phase of illumination from the 4 tetrahedron 
defining satellites, as shown by the thin dash-dot lines. 


the density of GaAs, the mass is about 52,000 kg/km^. 
The International Space Station, for comparison, is about 
420,000 kg, or equivalent to about 10 km^ of transmit¬ 
ter. It is possible that the mass of the lens might be 
reduced further, since thinner solar cells are in active de¬ 
velopment, though none yet have comparable efficiency. 
Thin film cells made of CdTe can have active areas as 
thin as 1 micron[5D], and there is research into cells just 
a few atoms thick, using the same high optical absorption 
exhibited by graphene PT]. 


Laser 

VCSELs (Vertical Cavity Surface Emitting Lasers) are 
a likely candidate for the laser amplifier. They can be 
built with the required efficiency, but require a lens since 
their native beam width is too high. 


Aiming the active lens 

Each pixel of the lens must know its location relative to 
the reference source and desired beam direction, in order 
to compute the correct phase. The coordinate system is 
defined by four spacecraft flying in formation, arranged 
in a tetrahedron a few hundred km on a side. Within 
the tetrahedron, the transmitter disk rotates, keeping it 
roughly flat. Each of the vertex spacecraft illuminates 
the transmitter disk. Each point on the disk compares 
the phase of the illumination, and its modulation, to 
determine its position. It then applies the appropriate 
phase correction to the reference signal, resulting in a 



thrust required 



FIG. 5: Stationkeeping for formation flying. (A) is the bal¬ 
ance of forces needed to say in the same relative position when 
at an angle Q from the central element, and at distance Ai? 
radially. (B) shows how reflective controls with micro-mirrors 
might work. 


single coherent output beam. The phase correction ap¬ 
plied will vary throughout the mission, both because the 
focal length is changing, and since the goal is not a single 
sharp focus, but a uniform intensity across the sail. 

The vertex spacecraft will need station-keeping. In 
particular, the vertex generating the reference beam 
needs to be almost directly in line with with the de¬ 
sired output beam and the center of the rotating lens. 
This is needed to keep the phase updates down to a rea¬ 
sonable rate. The other vertex spacecraft must at least 
measure, and potentially control, their location relative 
to the other vertices to within a fraction of a wavelength. 

The steering of the vertex spacecraft can potentially be 
done with solar sailing forces, as discussed later for the 
disk, or with tiny thrusters as developed for the LISA 
(Laser Interferometer Space Antenna) [52] mission. The 
measurement of the relative positions might use a combi¬ 
nation of the millimeter wave distance measurement de¬ 
veloped for GRACE[23] and GR AIL [5^ and their follow- 
on missions [25]. and the extremely precise optical mea¬ 
surement technology developed for LISA. 

The lens is a thin membrane and has very little struc¬ 
tural stiffness, though it will probably spin slowly to help 
keep it flat. Therefore we model it as an independent 
group of membrane elements that must fly in formation. 
In addition, the beam source (behind the lens) and the 
metrology elements (beyond the lens) must also fly in for¬ 
mation with the lens itself. What are the forces needed 
to do this, and how can they be generated? 

The main force that must be countered is the Sun’s 
gravity. This vector always points towards the Sun, 
and formation flying requires adding additional forces to 
make the vector force on each element equal to the force 
on the central element of the array. 

Gonsider an element that is at an angle B to the central 
element, as seen from the star. Let the distance from the 
central element to the star be i?, and the additional dis¬ 
tance for the element under consideration by AA. Now 
from the balance of forces as shown in Fig [^ we need 
a sideways of GM/R^ ■ sin(0), and an additional radial 
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Assuming at ±12° deflection from a static angle of 45°, 
then covering only 10“^ of the sunward surface with mi¬ 
cromirror arrays will generate sufficient thrust for sta¬ 
tionkeeping. Since these actuators require no power to 
remain in position, and cover only a small fraction of the 
surface, their effect on the mass and power budget should 
be small. 


Worked Examples 



FIG. 6: Micro-mirror array from Texas Instruments, from 
|26 | . (A) shows the operation of two pixels. (B) shows a por¬ 
tion of the roughly 1000 x 1000 array, with the pixels deflected 
one way in one image and the other way in the second image. 
Modern devices feature a ±12° degree deflection. 


force of GM{1/R^ ■ cos{9) — l/(i? ± Ai?)^. 

To get the required 4 kW on the array, it must be 0.58 
AU from the sun. If we wish to control arrays of up to 
100 km in radius, then 6* < 1.15 • 10“® radians. Assuming 
a maximum 45° tilt from the normal to the sun, then 
AR < 70.7 km. This implies a radial acceleration needed 
of 2.82-10“® m/sec^, and a perpendicular acceleration of 
2 • 10“® m/sec^. 

Consider now a 1 meter square of material. At 10 mi¬ 
crons in thickness, and a density of 5200 kg/m^, there is 
about 52 g of material. The incident solar energy is about 
4 kw, giving a force of 4000/c = 1.33T0“® Newtons. If we 
can convert this solar radiation force into steering thrust, 
we have plenty of authority for station keeping. We can 
do this by covering a small portion of the sunward side 
of the lens with steerable reflector arrays, such as the TI 
micro-mirrors pS] shown in figure 

One possible configuration could feature pairs of mir¬ 
rors, as shown in Fig 3(b). The normal un-deflected posi- 
ton would result in reflections b and e. Switching to c 
and d would increase the radial thrust, while switching 
to a and f would decrease it. Similarly, switching to a 
and d would generate a rightward thrust, while leaving 
the radial thrust largely unchanged. Reflections c and 
f would generate a leftward component. Other config¬ 
urations will generate a mixture of radial and sideways 
forces. Since each chip consists of millions of mirrors, and 
each mirror can be switched with microsecond speed, fine 
control of forces should be straightforward. 


Performance of the system is very dependent on the 
size of the transmitter. Increasing the radius r of the 
transmitter helps in two ways. As long as the tempera¬ 
ture limit is not reached, then acceleration scales as r^, 
plus the distance to which the beam can focus scales as 
r. Since, classically, v = \/2as, the final speed obtained 
scales as So a little less than 5-fold increase in ra¬ 
dius can increase the final speed by a factor of 10. We 
will therefore consider a “small” transmitter with a ra¬ 
dius of 7.65 km, a “medium” system with a radius of 32 
km, and a “large” system with a radius of 80 km. Assum¬ 
ing the mass of the transmitter is dominated by the mass 
of existing GaAs solar cells, this results in 25, 400, and 
2550 times the mass of the International Space Station 
(ISS). For comparison, the passive 1000 km diameter sail 
of Forward was 560,000 tons, or 1333 times the mass of 
ISS. 

The exact drive wavelength is not critical, though it 
should be far enough into the visible so that even when 
redshifted, it still remains in the region of high absorption 
as shown in Fig This requires the arriving photon to 
have an energy of at least 0.7 ev, and if the maximum 
speed obtained is 0.7c, then the source photon should 
have a wavelength of 745 nm or less. Here we assume a 
driving wavelength of 500 nm. 

First consider the small transmitter, pushing a 1 kg 
payload. This beam drives a 87.5m radius sail (also about 
1 kg), composed of 54 layer graphene. The sail absorbs 
72% of the incident light, reaching a temperature of 2865° 
K and generating 438 N of thrust. 

This thrust is conveyed to the 1kg payload by car¬ 
bon nanotube rigging (d). The thrust acting on the 
payload, sail, and rigging gives an acceleration of 219 
m/s^, or 22 gs. A 7.65 km aperture at 500 nm can fo¬ 
cus down to a 87m radius out to a distance of 14.7 AU. 
This allows the probe to accelerate for 1.6 days, reaching 
a speed of 0.103c. Relativistic corrections are small for 
such “slow” velocities, and including them gives a final 
speed of 0.100c. 

A bigger transmitter allows larger payloads, higher 
speeds, or both. Consider scaling the transmitter up to 
the medium size with a 32 km radius, still 3 times less 
massive than the passive lens in [T]. The first thought 
might be to drive the same 1 kg payload and sail with 
this larger transmitter, but that melts the sail. To keep 
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Example 

1 kg 

1000 kg 

4 kg 

1000 kg 

limited 

to 0.1c 

to 0.15c 

to 0.5c 

to 0.5c 

layers 

Radius of transmitter (km) 

7.65 

32.00 

32.00 

80.00 

32.00 

Area (square km) 

183.9 

3217.0 

3217.0 

20106.2 

3217.0 

Transmitter flux (w/m^) 

1000.0 

1000.0 

1000.0 

1000.0 

1000.0 

Power (GW) 

183.9 

3217.0 

3217.0 

20106.2 

3217.0 

Mass of transmitter (tonnes) 

9781.0 

171143.9 

171143.9 

1069649.5 

171143.9 

Mass in units of ISS 

23.3 

408.0 

408.0 

2550.1 

408.0 

Radius of sail (meters) 

87.5 

2767.8 

175.3 

2770.2 

7190.4 

Sail area (km^) 

0.024 

24.066 

0.097 

24.108 

162.426 

Number of layers 

54.0 

54.0 

54.0 

54.0 

8.0 

Sail mass (kg) 

1.0 

1000.7 

4.0 

1002.4 

1000.5 

Absorbence (%) 

71.54 

71.54 

71.54 

71.54 

16.99 

Sail absorb kW/m^ 

5466.3 

95.6 

23839.4 

596.6 

3.4 

Sail temperature (Kelvin) 

2865.1 

1042.0 

4140.4 

1646.8 

646.5 

Thrust (newtons) 

438.401 

7670.948 

7670.948 

47943.426 

1821.354 

Payload mass (kg) 

1.000 

1000.000 

4.000 

1000.000 

1000.000 

Acceleration, sail-|-payload only(m/s^) 

219.2 

3.8 

957.2 

23.9 

0.9 

Rigging mass (kg) 

0.000 

0.276 

0.017 

1.724 

0.170 

Acceleration including rigging(m/s^) 

219.1 

3.8 

955.1 

23.9 

0.9 

Wavelength of light used (nm) 

500.0 

500.0 

500.0 

500.0 

500.0 

Distance to cutoff (AU) 

14.7 

1941.1 

122.9 

4857.0 

5042.8 

V at cutoff, Newtonian (km/sec) 

31013.7 

47185.8 

187436.1 

186450.7 

37061.3 

Newtonian V (in term of C) 

0.103 

0.157 

0.625 

0.622 

0.124 

Newtonian time (days) 

1.6 

142.5 

2.3 

90.2 

471.2 

Relativistic V (in term of C) 

0.100 

0.149 

0.502 

0.500 

0.119 

Relativistic time (days) 

1.7 

146.3 

2.5 

100.6 

481.1 


TABLE III: Worked examples of different performance levels 


the temperature down, we need at least a 4x bigger sail 
(and even that gives 4140° K). This accelerates the 4 kg 
payload at a little less than 100 g. A classical calcula¬ 
tion gives a final speed of 0.625c and hence is clearly not 
appropriate. The relativistic calculation of section III 
shows a final speed of 0.5c. 

Alternatively, also with the medium size transmitter, 
we can launch a larger payload more slowly. A 1000 kg 
payload, with a correspondingly larger sail, accelerates 
at only 0.4 g. However, we can keep this up for 146 
days, leading to a final velocity of 0.149c after relativistic 
correction. The sail in this case is cool at less than 1100° 
K. 

The fourth example demonstrates that we can use an 
even larger transmitter to send big payloads fast. An 80 
km radius transmitter can send a 1000 kg payload to 50% 
of the speed of light in 100 days at 2.5 g. 

The final example shows that even if the extrapolated 
absorption trend does not continue to 54 layer graphene, 
all is not lost. It uses the same 32 km transmitter and 
1000 kg payload as the second example, but limits the 
number of layers to the experimentally measured range. 
The final attainable velocity is decreased only slightly. 


from 0.149c to 0.119c. A bigger drawback is that the 
roughly 3-fold lower acceleration means the time to final 
velocity is that much greater. 

GETTING THE DATA BAGK 

Assuming we can launch a flyby probe, and it collects 
data as it passes through an alien solar system, we then 
have the problem of how to get the data back. Solar 
energy from the target star will be insufficient, as the 
probe is only within any substantial illumination for at 
most a few hours. 

The large transmitter provides a potential answer to 
this dilemma. The gain of a uniformly illuminated, co¬ 
herent emission disk is 

G = (27rr/A)2 

Multiplied by the large transmitted power, this results 
in a huge EIRP (roughly 5.20 • 10^^ watts for the 32 km 
radius transmitter, compared to the roughly 3.85 • 10^® 
watts of the sun). This is enough to provide 4.6 w/m^ 
at a distance of 10 light years. Since this will arrive 
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in monochromatic (though redshifted) form it should be 
especially efficient to convert this to useable electrical 
energy. 

Now to send the data back, assume a 1 meter aper¬ 
ture, and a 1 watt optical transmitter at A = 500 nm. 
Then the EIRP, computed as above, is 3.95 • 10^^ watts. 
After 10 light years, this gives a received power density 
of 3.51 • 10“^^ w/m^. Since each photon has 3.98 • 10“^® 
joules, a large collector is clearly needed to collect a sig¬ 
nificant number of photons. Since we are assuming large, 
space-based structures are possible, perhaps assume a 
receiver about a quarter of the area of the transmit¬ 
ter. Such a 16 km radius telescope would collect about 
700,000 photons per second. If we can pick them out of 
the background light from the star, this should give a de¬ 
cent download rate. The star (if sunlike) is roughly 10^^ 
times brighter. However we can apply very aggressive 
optical filtering (since the data light is monochromatic 
and of known wavelength), temporal filtering (we could 
perhaps use 1 ps pulses every microsecond, for another 
gain of 10®), and coronograph techniques currently being 
developed for the search for extrasolar planets. Between 
these we should be able to distinguish the data photons 
from the noise and get the probe’s data back. 

OTHER USES 

Flyby missions for an Oort cloud object would be easy 
and cheap (no on-board power needed). 

Such a system could make short work of any orbital de¬ 
bris, and defiect incoming asteroids or cometsP?]. Con¬ 
versely, it would make a fine weapon system, dumping 
enormous power into a tiny area on any solar system 
body. This possibility might argue for putting the trans¬ 
mitter outside the orbit of the Earth. This would require 
a larger transmitter for the same output, but be physi¬ 
cally unable to target the Earth. 

Such a system could make a completely unambiguous 
way to make contact with an extraterrestial civilization, 
provided only the arriving wavelength is within the vis¬ 
ible spectrum of the observer. An observer in the beam 
would see the parent star suddenly increase in brightness 
by 10® times or so, until (for nearby stars) it casts an ap¬ 
preciable fraction of daylight (at 10 light years, it would 
be considerably brighter than typical office lighting here 
on Earth). Any society that watches the sky could not 
miss such a signal. A message could be sent by blink¬ 
ing such a beam at a slow rate (perhaps 1 bit/sec) in a 
pattern that repeats with a period that is a product of 
two primes[5S]. Any society that has discovered prime 
numbers and factoring (such as the ancient Greeks or 
Chinese) might note that the repeat interval is the prod¬ 
uct of primes, say M x N, and then display the data on a 
M X N grid, and then attempt to decipher the message. 
Such a system could therefore contact a civilization that 


is intelligent but not technological. 

FUTURE WORK 

The required power, and transmitter size, could in the¬ 
ory be further reduced by propelling the sail with UV 
light, where the absorption of a single layer of graphene 
reaches 10% [2^. However, there is currently no known 
method for generating UV light with high efficiency. 

DISCUSSION 

Graphene sails and phased array optical drive can po¬ 
tentially reduce the scale of structures needed for inter¬ 
stellar probes. No new physics is needed - “only’ the 
technologies of multi-acre graphene sheets, and kilometer 
scale nanotube ropes, and the ability to build large space- 
based structures. For comparison, the LHG is about 4.3 
km in radius, the ISS about 420,000 kg, and large photo¬ 
voltaic power plants are perhaps 5 km^ in size and more 
than 0.5 GW in output. These are within 1-2 orders of 
magnitude of what is needed. Furthermore, most of the 
technologies needed (efficient solar cells, efficient lasers, 
precision space constellations, graphene and nanotube 
production) are being developed for other reasons. 

For such a system, the vast majority of the cost is in 
the transmitter. The probes are themselves should be 
relatively cheap. Furthermore, the launch time is rela¬ 
tively short, then there is a long coast where the time 
depends on the destination, then a relatively quick data 
retrieval. So this architecture makes sense for launching, 
then retrieving data from a fleet of probes. 

CONCLUSION 

This paper continues the work of investigating different 
possible architectures for interstellar probes, searching 
for those that may be more economical and practical. 

Even with the possible improvements outlined here, 
the facilites for interstellar probes will still be large, and 
too expensive for the current earth economy. They will 
not be built in the next few years, nor in the next few 
decades, but perhaps they are not so far out of reach as 
has been thought. 

Appendix 1: Code for calculations 

Here is a program in C-I--I- that performs the calcula¬ 
tions used in the table. The input is a file that contains 
examples, one per line, with the transmitter radius in me¬ 
ters, the maximum number of layers in the sail, the areal 
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power output of the transmitter, and the payload mass. 
The output is a table in LaTex format, called “tab.tex”. 

This program does not use the analytic results from the 
paper concerning the number of layers or the optimum 
payload/sail mass ratio. Instead the optimum is found 
by gradient descent. 

#include <stdio.h> 

9 ^include <math . h> 

#include <vector> 
using namespace std; 


const 

double 

GRAPHENE^BSORBS = 

0.0230; 

// 

one layer 


absorbs 

2.3% of visible lig 

h t 



const 

double 

GRAPHENE-MASS == 

0.77e-6; 

// 

mass of 0.77 mg 


/m-2 





const 

double 

GRAPHENE-THICK = 

0.345e-9; 

// 

thickness of a 


layer 





const 

double 

GRAPHENE_STRENGTH= 

1.0 el 1 ; 

// 

Pascals , 


graphen. 





const 

double 

RIGGING-STRENGTH = 

100e9 ; 

// 

Pascals , carbon 


n an 0 t u 

be 




const 

double 

RIGGING-DENSITY = 

1300; 

// 

kg/m~ 3, carbon 


n an 0 tub 

e s 




const 

double 

SI-DENSITY = 

2329; 

// 

k g /m ~ 3 , silicon 

const 

double 

GaAs-DENSITY = 

5320; 

// 

k g /m ' 2 , Gallium 


Arseni. 

de 




const 

double 

STEFAN-BOLTZMANN = 

5.67e-8; 

// 

Stefan- 


Boltzmann constant 




const 

double 

C = 

3e8 ; 

// 

speed of light , 


m/s 





const 

double 

G = 

9.8; 

// 

g at Earth 


surface 

m/s ~2 




const 

double 

AU = 

1.496 el 1 ; 

// 

AU in meters 

const 

double 

LAMBDA 

= 500e-9; 

// 

500 nm light 

const 

double 

LENS-THICKNESS 

= lOe-6; 

// 

10 microns 


thick 





const 

double 

ISS-MASS 

; 419455; 

// 

mass of 


interna 

tional space station 

in kg 



class 

print 

{ 





public ; 

const char >i>label ; 
const char ^format ; 
double scale ; 
bool line ; 

print(const char *a, const char >kb = double s = 

1.0, double 1 = false){label = a; format = b; scale 

= s ; line = 1 ;} 

}; 


// calculate one specific case 

double one_calc (double tr , double sr , int nl , double pmass , 
double p o w e r _p e r _ a r e a , v e c t o r <doub le > &rslt ) { 
rslt[0] = tr; 
rslt [6] = sr; 
rslt [8] = nl; 
rslt [14] = pmass; 

printf{”\nTrans_radius_%.lf ,_sail_rad_%.lf , _%d _ 1 a y e r s , _ 
p ay 1 o a d _ % . 3 f \ n ” , tr , sr , nl , pmass) ; 
double tr_area rslt [1] = M_PI *pow(tr ,2) ; 

rslt [2] = power_per_area; 

double power = rslt [3] = tr_area * power_per_area; 

printf (•• area_of_transmitter_%.l f _km ~ 2 ,_power_%.lf_gigawatts\ 
n”, tr_area/le6, power/le9); 

double tr-mass = rslt [4] = rslt [5] = tr.area* LENS-THICKNESS 

*GaAs_DENSITY ; 

printf (”Mass_of_lens_%.lf_kg,_%.2f_space_stations\n” , 
tr_mass , t r _ m a s s / ISS _M AS S ) ; 

double absorb = rslt [10] = 1 - pow(l-GRAPHENE-ABSORBS, nl); 
printf{”A_sail_ wi th_%d_1 ayers _absorbs . 2 f%%_ of_incident _ 

light.\n”, nl, absorb*100); 
double sail_area = rslt [7] = M_PI x- pow(sr ,2) ; 

double s a i 1 _ p o w e r _ p e r _ a r e a = rslt [11] = power * absorb / 

sail_area ; 

double sail-mass = rslt [9] = sail-area * GRAPHENE-MASS * nl; 

double temp = rslt [12] = pow ( sail-power-per-area /(2* 

STEF AN-BOLTZMANN* absorb ) , 0.25) ; 

printf(” Sail_area_%.lf _km '2,_power_density_%.lf _w/m' 2 , _temp _ 
% . 1 f _K , _mass _% . 3 f _kg\n” , 

sail-area/le6, sail-power-per-area, temp, sail-mass); 
double thrust = rslt [13] = s a i 1 - p o w e r - p e r - a r e a * sail-area / 

C; 

double acc = rslt [15] = thrust / (pmass + sail_mass); 

printf(” Initial_guess_at_acceleration_%.3f _m/ s,_%.3f_gs\n” , 
acc , acc /G) ; 

double r i g g i n g-1 e n g t h = sqrt(2) * sr; 

double r i g g i n g - c r o s s - a r e a = sqrt(2) * pmass*acc / 
RIGGING-STRENGTH; 

double rigging-mass = rslt [16] = RIGGING-DENSITY * 

rigging-length * rigging-cross-area ; 
p r i nt f ( ” Mass_of _ r i gg i n g _%.3 f _kg\n” , r i g g i n g - m a s s ) ; 
acc = rslt [17] = thrust / (pmass + sail-mass + rigging-mass); 

printf (” Better_guess_at_acceleration _9o.3f _m/ s,_%.3f_gs\n” , 
acc , acc /G) ; 
r s 11 [ 1 8 ] = LAMBDA; 

double dist = rslt [19] = 4 * sr * t r/LAMBDA/2.4 4; 

printf (” Distance_until_can ’t_focus_9o.lf _km . . 2 f _AU\ n” , dist 

, dist/AU) ; 

double V = rslt [20] = rslt [21] = sqrt(2 * acc * dist); 
double t = rslt [22] = sqrt(2 * dist / acc); 

printf (” Newtonian_v = %.lf _m/ s , _%.4f_c ;_t=%.lf_sec , _%.2 f_days 

_%.3f_years\n”, v, v/C, t, t/86400, t/365.25/86400); 
double err = 10; 


do { 

double eterm = —expml(—acc/C^t) ; 

V = C * eterm ; 

double s = C* t — C*C/acc * e t er m ; 
err dist—s; 

t +— err/v; 

pr i ntf ( ” Re 1 at i vist i c _s=%.l f _v = %,l f _er r =%.3f\n” , s, v, 
err) ; 

} 

while (fabs(err) > 100.0); 

rslt [23] = v; 

rslt[24j = t; 

printf(” Relativistic _v = %.l f _m/s , _%,4 f _c ;_t=%.l f _sec , _%,2 f _ 

days_%.3f-years\n”, v, v/C, t, t/86400, t/365.25/86400) 

double thick = n 1 *GRAPHENE-THICK; 
double hoop = GRAPHENE-STRENGTH * 6.0/13.0; 

double edge = rslt [25] = t h r u s t / ( 2 >i< M-PI ) / t h i c k / hoop ; 

printf (” Requires_outer_%.3f_meters_to — make _hoop — stress — 
positive\n”, edge); 
return v/C; 

} 

int main(int argc , char X'X'argv) { 

printf(”enter_radius—of — transmitter_(meters) ,—radius—of — sail 
— (meters) ,—number — of — monolayers , —mass — of — payload — (kg) :\ 
n” ) ; 

V e c t o r <doub le > trad , srad , pm, ppa ; 
vector<int> layers ; 
for(;j) { 

double tr , sr , pmass, power-per-area; 

if (scanf(”%lf —%d—% 1 f _% If” , &tr , &nlmax, &pmass , & 

p o w er_per_area ) != 4) 

break ; 

printf ( ” Got-%f\n” , tr); 
trad.push-back(tr); 

srad.push-back(tr/10.0); // start with sail 1/10 

transmitter size 

layers.push-back(nlmax); // max number of layers 

pm . push-back (pmass) ; 
ppa.push-back(power-per-area); 

} 

ve c t or < V e c t o r <double> > table ; 
for(int i=0; i <t r ad . s i z e ( ) ; i++) { 

int LMAX = 1 a y e r s [ i ] ; 

1 ay e r s [ i ] = 1 ; 

/ / loop until find the right value 
ve c t o r <double> numbers(26) ; 
for(;;) { 

V e c t o r <double> rslt ; 

double delta = 0.001 * srad[i]; 

for(int n^layers [ i ] — 1; n <= layers[i] + l; n++) { 

for (double r = srad[i] — delta; r <srad[i] + 1.5* 

delta; r += delta ) 

rslt.push_back(one_calc(trad[i], r, n, pm[i], 
ppa[ i ] , numbers) ) ; 

} 

if (layers [i] + l > LMAX) 

rslt [6] = rslt [7] = rslt[8] = 0.0; 

pr i nt f ( ”_rad=%9.1 f-%9.1 f-%9.1 f \n” , srad[i]-delta, 

srad[i] , srad[i]+delta); 
for(int n = 0; n< rslt . size () ; n +=3) 

printf (”n=%3d-%9.5 f-%9.5f-%9.5 f \n” , layers [ i ] + 

(n/3-l), rslt[„], rBlt[n + l], rslt[n + 2])j 
//find the biggest one . If it’s not numb er 5 , 
readjust and try again. 
double biggest = 0.0; 
int bigi = —1; 

for(int k=:0; k<rslt.size(); k++) { 

if (rslt [k] > biggest) { 

biggest = rslt[k]; 
bigi = k; 

} 

} 

printf (” Biggest — is—case —%d\n” , bigi); 
if (bigi == 4) 
break; 

int ncase = bigi/3; 

layers [ i ] += (ncase —1) ; 

int rcase = (bigi % 3); 

srad [ i ] += (rcase — 1) >k delta; 

} 

// OK, now (re) compute final result 

one-calc(trad[i], srad[i], layers[i], pm [ i ] , ppa [ i ] , 
numbers ) ; 

table . push_back (numbers) ; 

} 

print labels [] = { p r i n t ( ” Rad i u s — o f - t r a n s m i 11 e r — ( km ) ” , ” 9o . 2 f ” 

. le-3), // [0] 

print (” Area—(square —km )” , ” % . 1 f” , 1e 

- 6 ), // 111 

print (” Transmitter — flux — ( w/m$ ' 2 $ ) ” , ” 

%.lf” , 1.0) , 

print ( ” Power - (GW) ”, ” %.l f ” , le - 9) , 

// [SI 

print (’’Mass — of — transmitter—(tonnes)” , 
”%.lf” , le-3) , 

// [ 6 ] 

p r i n t ( ” Mass - i n - u n i t s — o f — IS S ” , ”%.lf” , 

1,0/ISS-MASS , true), // [7] 

print (” Radius-of—sail —(meters)” ) , 

// ISI 

pr int ( ” S ail „area_ (km$'2$ ) ” , ”%.3f”, le 

- 6 ), // [91 

print(”Number-of—layers”), // 

[ 4 ] 

print (” Sail-mass-(kg)”) , // [10] 

print(”Absorbence —(\\%)” , ”%.2f” , 

100.0), // [8] 
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print {” S ail _absorb_kW/m$''2$” , , 

le-3), // [11] 

print (”Sail_temperature_{Ke Ivin)”), 

// [12] 

print {”Thrust_(newtons)” , 1,0, 

true) , // [13] 

print (” Payload_mass_(kg)” , ”%.3f” , 

1.0) , // [5] 

print (” Acceleration , _sail+payload_only 

(m/s$'2$)”) , // [14] 

print (” Rigging_mass_(kg)” , ”%.3f”) , 

/// [15] 

print (” Acceleration_including_rigging ( 
m/s$''2$)”) , // [16] 

print ( ” W avelength_of_liglit _used_ (nm) ” , 
, le9) , // [1 7] 

print {” Distance_to_cutoff_(AU)” , ”%.lf 

" , l.O/AU) , // IIS] 

print (”V_at_cutoff ,_Newtonian_ (km/ sec) 
”, , le-3) , 

// 19 

print (” Newtonian_V_ ( in _term_of_C) ” , ” 

%.3f”, 1.0/C), // 20 

print ( ” Newtonian _t i me _ ( d ays ) ” , 

1.0/86400.0), // 21 

print (” Relativistic _V_ { in_term_of _C) ” , 
”%.3f” , 1.0/C) , // 22 

print {”Relativistic_time_(days)” , ”%.l 

f”, 1.0/86400.0), // 23 

print {” Rim_width_ ( meters )” , ” %.3 f ” ) } ; 

// 

// write table in tex format 
FILE »ftab = fopen(”tab.tex”,” w” ) ; 
if ( ft ab === NULL) { 

printf(”can’t_open_tex_file\n”); 
return 4 2; 

} 

const char *cs = (argc > 1) ? ”%” : ””j y / comment string 


fpri 

nt f ( ftab 

, ”%s \\ document cJ 

.ass[10pt]{ article }\n” 



”%s\\usepackage 

{array}\n” 



”%s\\usepackage 

[margin=:0.5in]{geometry}\ 



”%s\\begin{docu 

ment}\n” , cs ,cs ,cs ,cs) ; 

fpri 

nt f { ftab 

, ”%s\\begin{tabu 

ilar}{_r” , Cf.) ; 

for { 

int i=0; 

i<table.size(); 

i++) 


fprintf (ftab , ”_r”) ; 


fpri 

nt f ( ftab 

, ”}\n”); 


fpri 

nt f ( ftab 

, ”%sExample_ID_” 

, c=) ; 

for { 

int i=0; 

i<t ab1e . size () ; 

i++) 


fprintf(ftab , ” &c- (% c ) ” , 

■A' + i); 

fpri 

nt f { ftab 

, ”\\\\_\\hline\n 

”) ; 

for ( 

int i=0; 

+ 

+ 

II 

V 



fprintf(ftab , ”%s_”, labe 

lE, [ i ] . labal ) ; 


for(int 

k:=:0; k<table . size 

() ; k++) { 


fpri 

nt f ( ftab , ) ; 



fpri 

ntf ( ftab , labels [ 

i ]. format , t a b 1 e [ k ] [ i ] » 1 a 



].scale ); 



fprintf(ftab , ”\\\\%s\n” , 

labels [i], line ? "\\hlin 


”) ; 


} 


fprintf (ftab , 

” \\end{t abul 

ar}\n”); 


fprintf (ftab , 

”%s \\end{doc 

ument }\n” 

, cs 


£c los e ( ftab ) ; 
return 0; 

} 
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